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Abstract 

When a single beam-splitter receives two beams of bosons described by Fock states (Bose-Einstein 
condensates at very low temperatures), interesting generalizations of the two- photon Hong-Ou-Mandel 
effect take place for larger number of particles. The distributions of particles at two detectors behind 
the beam splitter can be understood as resulting from the combination of two effects, the spontaneous 
phase appearing during quantum measurement, and the quantum angle. The latter introduces quan- 
tum "population oscillations", which can be seen as a generalized Hong-Ou-Mandel effect, although 
they do not always correspond to even-odd oscillations. 
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Beam splitters are an essential component of many experiments designed to observe quantum effects. 
They are involved in experimental and theoretical schemes that both Helmut Ranch and Daniel Green- 
berger have studied. Indeed, the famous neutron experiments of H. Rauch and colleagues [11 [2] were made 
possible by the realization of an appropriate device allowing neutron beams to be split into two coherent 
beams, which can then be recombined and give rise to various interesting quantum interference effects. 
The observation of equally famous quantum GHZ (Greenberger, Horne and Zeilinger) violations of local 
realism [31 0] may also require the use of photon beam splitters ^ . Still another example is given by the 
entanglement swapping effect, which requires indistinguishable photons to be measured at the output 
ports of a beam splitter [6 . This list is, of course, non-exhaustive. 

Here, we come back to the basic properties of a single beam splitter and show that, as simple as it 
may look, it already exhibits strong quantum properties. Previous studies of the quantum properties of 
beam splitters include Refs. [TUSHH]- Generally, beam splitters are used in conditions where they receive 
particles one by one. Here we generalize the discussion and consider the case where a beam splitter 
receives groups of particles in its two input beams, described by Fock states of bosons. We will study 
the effects of the "quantum angle", which was introduced in the context of more elaborate interferometry 
experiments involving several beam splitters, and leads to violations of local realist BCHSH and GHZ 
inequalities |10l [TT] . Holland and Burnett have studied the quantum limits on the detection of small 
phase shifts with interferometers involving two beam splitters [12]; for this purpose they also study the 
distribution of the relative phase of the two output beams at a single beam splitter, assuming that the 
two incoming Fock states have equal populations (twin states) . Here we release this assumption and study 
the distributions of the immber of particles at the outputs. 

The production of Fock states with photons is not an easy task, if not impossible, except for a 
small immber of photons; see for instance Ref. |13| for a description of an experiment with two states 
containing two photons. Coherent states are, of course, much easier to produce, even with a small 
average number of photons (very small intensities), but they remain fundamentally very different from 
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Fock states. Fortunately, the phenomenon of Bose-Einstein condensation in ultra-cold gases provides us 
with a method to produce condensates thermally and, when repulsive interactions between the atoms 
stabilize the condensate, there are good reasons to believe that its state is well described by a Fock state. 
Thermal excitations are of course always present, but they can be reduced very efficiently by reducing 
the temperature. Moreover, the technique of Bragg scattering of atoms from standing laser waves [Hj can 
be used to obtain efficient atomic beam splitters |15j . and even observe interferences with Bose-Einstein 
condensates in interferometers with macroscopic arm separation |16| . 

1 Classical and quantum calculation 

The situation we consider is shown schematically in Fig. [1] We will perform a quantum calculation but, 
as a point of comparison, we start with a simple classical calculation. 
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Figure 1: Na^Np bosons proceed from the sources to a beam splitter, followed by two detectors 1 and 
2, where mi and m2 particles are detected. A phase shift is inserted in the a-arm for generality, but 
turns out to play no role in the results. 



1.1 Classical model 

In classical optics, if two beams with equal intensities Iq are sent to a beam splitter with a relative phase A, 
the output intensities Ii and I2 at the output detectors 1 and 2 are proportional to [1 + cos (A — 7r/2)] and 
[1 — cos (A — 7r/2)] (the 7r/2 arises because a phase shift occurs at a reflection, but not at a transmission). 
If the phase A is completely unknown, these expressions have to be summed over A between — tt and +7r; 
a well-known classical calculation then shows that the distribution P{I) of the random variables /1.2 is 
given by: 

If the input intensities are different, la and Ip — x^Ia, this calculation can easily be generalized. The 
two intensities are now proportional to [1 + r cos (A — 7r/2)] and [1 — r cos (A — 7r/2)], where: 
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These expressions result from purely classical wave theory. 

Semi-classical expressions can be obtained by considering a flux of classical particles reaching inde- 
pendently the beam splitter, each having a probability [1 + rcos (A — 7r/2)] /2 to go to detector 1, and a 
probability [1 — r cos (A — 7r/2)] /2 to go to detector 2. The probability that, among a total of N particles, 
mi will go to detector 1 and m2 to detector 2 (with mi -I- m2 = N) is then given by: 

TV' 

class, (mi, ma) = TTXTTT^TZ-r [1 + ^ cos (A - 7r/2)]™^ [1 - r cos (A - 7r/2)]"^ (4) 
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For Fock states, we expect that the relative phase A should be completely random, so that this expression 
becomes: 

Xmi,m2) = T— 7 / — [l + rcosA]"Ml-rcosA]"' (5) 
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Figure [21-a shows an example of such a distribution for equal intensities of the incoming beams, which 
reproduces the shape of the classical distribution ([1]), with a minimum at the center and maxima at the 
edges. Figure [2]-b shows another example, assuming now that the intensities of the two incoming beams 
are different, and that their ratio is 6/44. Because the interference effect between the two beams can no 
longer be completely destructive, the distribution tends to concentrate more towards the center or the 
curve (medium values of mi). 
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Figure 2: The left part (a) shows the classical (dotted line) and semiclassical (full line) distributions as 
a function of mi when the intensities of the input beams are equal {x — r ~ 1) and when mi + m2 = 50. 
Values of mi near the maximum and the minimum are more likely to occur with this distribution. The 
right part (b) shows the same distributions for the same total number of particles, but when the intensities 
of the input beams are different (their ratio is 6/44). 



An equivalent summation formulation more suitable for computations is found by expanding the 
binomials (1 ± rcos A)™' and integrating term by term. The result is: 
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where K is a. normalization factor. This formula was used, for instance, to produce Fig. [5] and the dotted 
hne of Fig. Hb. 
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1.2 Quantum calculation 

We give two calculations, for input beams described either by Fock states or by coherent states. 



1.2.1 Fock states 

Before the beams of bosons cross the beam splitter, they are described by the quantum state: 

\N^,N0)^^^==ai^-ap\O) (7) 

Our calculation that is essentially the same as that of [lOj and [llj . The destruction operators associated 
with the two output beams (and detectors) are: 

ai = {e'-^Qa + iap) ; 02 = [ie^'^aa + ap) (8) 

The amplitude for finding mi,m2 particles in the detectors given sources with Na,Np particles in the 
sources is: 

^ ^NJNpl e'x mil ma! 

V^^m!^ (72)"'+™—* p!(mi-p)!(z!(m2-g)! 

X Sp+q,Na,^mi+m2-p-q,Nff (9) 

where x is a phase factor without physical relevance. Two methods of calculation are now possible. 
It is possible to replace the second (5-function in Eq. ^ by: 

X M - [ ^^^p#("ll+™2-p-<2-A'^^) (ir\\ 

to obtain: 

X (e^«+je'^)'"^(ze"' + e^*)"^ (11) 

The square the modulus of this expression contains an integral over two variables ip and (p'; if we make 
the changes of variables: 

; A=*^ (12) 

we find for the probability the expression: 

P(mi,m2) = ^^ — — cos [{N^-N^) A] 

[cos A + cos A] [cos A - cos A] (13) 

(note that the phase shift 9 has disappeared from this result). Assume for a moment that the A can be 
replaced by in the three cosines that contain it. Then A disappears, and we are left with an expression 
that is identical to (O with r = 1, except for normalization factors. We therefore see that A (or, more 



4 



precisely, A + 7r/2) plays the role of the classical relative phase of the two sources; this phase is averaged 
over 27r, which is normal since the phase in a Fock state is completely undetermined. For this reason, we 
will call A the classical phase angle, and A the quantum angle; we study in more detail below how the 
non-zero values of A introduce quantum effects. 

Another method is to use the 5— functions to eliminate the summation variable q in Eq. ^ and then 
square the result. We then find: 
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This expression is more convenient than (1131) for accurate numerical calculations. 
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1.2.2 Coherent states 

We replace the ket ([7]) by a product of coherent input states: 
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where E^ and Ef^ are complex number defining the intensity la.p — \Ea^f}\^ and the phases ^Pa.p of the 
incoming fields. In this expression, we can replace the creation operators by their expressions obtained 
from ([5]) and obtain: 
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Therefore, the output state is a product of coherent states as well, with amplitudes of the fields given by 
{Ea + iEp) /\/2 and [lEa + Efi) /\/2, which correspond exactly to the classical formulas. The operators 
contained in the exponential commute. By expanding it into a series as in (1151) . we obtain the probability 
to measure mi bosons at output 1 and m2 at output 2 as a product of Poissonian distributions: 
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[1 + r cos (A - 7r/2)]"^ [1 - r cos (A - tt/2)X 



where A = (/Jq — + 0. The result is therefore very similar to Eq. as well as to ([5]) if we assume that 
the phases of the incoming coherent beams are random. 

1.3 The generalized beam splitter theorem 

The properties of the distributions obtained in (ITSl) and yield the generalized Hong-Ou-Mandel 
theorem as we see next. 



1.3.1 Calculation 

The second line in Eq. (ITU) can be factored to produce: 
where: 



= ( cos- j ( sin- 



(18) 
(19) 
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with: 

= 0-0+1 (20) 

An exact result then is: 
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If Na = = 1, we are in the situation of the Hong-Ou-Mandel effect [T7|. When two photons fall on 
a beam splitter from two symmetrical directions, it is known that a quantum interference effect prevents 
them from leaving the beam splitter separately; they always leave together (in the same direction) . Here 
we obtain a direct generalization of this theorem: if an equal number of particles approaches from each 
side to meet at the beam splitter, an even number must emerge from each side. This is because, if 
Na = Np in the probability amplitude (j2ip . we have N = mi + m2 even, in which case mi and m2 are 
both even or both odd; but, if m2 is odd, the integral is over an odd function and therefore vanishes. 

1.3.2 A Gaussian fit 

A plot of the second line in Eq. (PT|) is shown in Fig. [31 We see from the figure that a (double) Gaussian 
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Figure 3: A plot of Q{(f)) for mi = 17 and m2 — 83 and the phase — tt/2. The peaks are at 00 = ±0.73 tt 
(the phase choice gives symmetrical peaks about zero). The relative sign of the two peaks is (~1)™2. 
The peaks are normalized to unit height here. 

fit to (5(0) is likely to be a good approximation. We set the first derivative of the logarithm of (5(0) to 
zero, which gives the maxima positions at ±0o where: 

0o = 2arccos(^y^^ (22) 

and also find that the second derivative there is —2N. The peak at — 0o is negative if m2 is odd. Hence 
we get the approximation: 

g(^) = QY''^ {^r'" [e-(^-^°)^ + (_l)™.e~f (^-^^)^] (23) 

where the prefactor to the Gaussians comes from Q{(f)o) upon use of cos^ 0o/2 = rrii/M and sin^ 0o/2 = 
m2/M. In Eq. (|2ip we then Fourier transform the Gaussians; doing these integrals and squaring gives 
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the probability: 
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for 7712 even 
for m2 odd 



An interesting further approximation uses the Stirhng formula for the to.;!. We have 



(24) 



(25) 

The factor of 1/2 in the exponent is usually dropped, but is actually important in our case in giving a 
characteristic shape to the probability curves. Moreover it makes the Stirling formula accurate to within 
a few percent for rrii > 2. The result is then: 
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Note the characteristic ^mim2 form in the denominator stemming here from the Stirling formula. This 
formula matches the exact result in Eq. (fT4)) very accurately for > 2. 



2 Physical discussion 
2.1 Few bosons 

We now study in more detail the consequences of the rule according to which even numbers of particles 
emerge from the beam splitter if an equal number impinge on each side. For this purpose, we use Eq. 

to calculate the distribution corresponding to the various possible numbers of particles detected at 
1 and 2. 

Prob 




Figure 4: The probability for Nq. = 4, Np — 4 illustrating the rule that, if an even number of particles 
enters each side of the beam splitter, an even number must emerge from each side. 
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Fig. m shows four particles entering each side of the beam spUtter. As expected, only an even number 
can emerge on each side, which explains the zeroes of the curve. Comparing with Fig. [5] immediately 
indicates that these cancellations superimpose strong quantum oscillations onto the classical intensity 
distribution. We have a sort of combination of a classical average over a phase A with rapid variations 
created by the Hong-Ou-Mandel effect. 




Figure 5: The probability for Na = 4, Np = 5 . The probability no longer vanishes for odd mi but 
oscillations remain. Now, of course, odd mi implies even m2 and vice versa. 

For Na — 4, — 5, the rule no longer applies, but the calculation of the distribution can still be 
done. The result is shown in Fig. [SJ which again contains odd-even variations and oscillations, even if 
the probability does not vanish for any value of the m's. In this case also, we have the superposition of 
an average curve, which can be understood in terms of a classical phase, and of an oscillation that can 
be seen as a generalized Hong-Ou-Mandel effect. 



2.2 More particles 

For larger values of N, the even rule is shown in Fig. [SI where the variation of the probability with mi is 
plotted for Na = Nf^ = 25. The characteristic probability variation l/^mim2 given in Eq. is visible, 
but again with strong quantum oscillations. In this case in Eq. (j26D the factor sin [(A^q — A^/3)0o/2]^ 
vanishes to satisfy the even rule. 

Suppose now we have slightly different source populations, Na — 26 and Np — 25. The result is shown 
in Fig. [71 When Na — Np = 1 we can find explicitly that: 
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The probability oscillates between these two values, modulated by the l/y^mimj factor, and changes 
over from maxima at even to odd at mi = 25. It is interesting to see that, in this case, the quantum 
oscillations vanish at the center of the distribution, but remain very pronounced at both sides. 

Larger population imbalances in the sources, at constant sum N = 50, result in even more complicated 
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Figure 6: The probability for Na — Nfs — 25 . The probabiUty vanishes for odd mi. The graph has 
been cut off at mi = 0, 50 where it is about twice as high. The variation with l/y/mim2 is evident. The 
dotted hne shows the corresponding semi-classical distribution 
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Figure 7: The probability for = 26, Np = 25 



behavior. What happens for instance if Na 



Nfs = 2 is shown in Fig. [51 We can then show that: 
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The probability oscillates between these two curves with "nodes" at mi = 7 and 43 corresponding to the 
crossing of the two quadratic curves. The nodes are actually a consequence of the discrete character of 
mi: if mi is replaced by a continuous variable in Eq. (I^Tl) . then the probability distribution becomes 
an oscillating function with a slowly varying amplitude. If the maxima and minima occur near integer 
values of mi, they remain very visible in the discrete version of the distribution, resulting in antinodes; 
but, if they occur near half integer values, the oscillations disappear in the discrete version, resulting in 
nodes. This "stroboscopic effect" also explains the minimum of oscillations at the center of Fig. [T] 

The "beating wavelength" becomes shorter as Na — Np becomes larger. The case Na = 28, Np = 22 
is shown in Fig. I^l-a; the case Na — 44, Np = 6 is shown in Fig. [S]-b. In this case, the fast oscillations 
of the generalized Hong-Ou-Mandel type have disappeared and have become significantly slower, which 
presumably makes them easier to observe experimentally. 
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Figure 8: The probability for TV^ = 26, Np = 24 
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Figure 9: The left part (a) shows the probability distribution for Na — 28, Np — 22. The right part (b) 
shows the same curve for Na = 44, Np — 6 (full line) and, as a point or comparison (dotted line) the 
classical curve of Fig. ^h. 



2.3 Role of the quantum angle 

If Na — Np, we can compare the semi-classical expression, Eq. ([S]) with r = 1, to the quantum result 
Eq. (fT5|) . We see that, instead of classical probabilities [IzbcosA], the quantum expression contains 
quasi-probabilities [cos A ± cos A] , which can take negative values when the quantum angle A does not 
vanish. This introduces quantum effects in a way that is reminiscent of quantum effects arising from 
negative values of the Wigner transform. The quantum angle is responsible for the population oscillations 
introduced by the beam splitter. 

When Na 7^ Np, Eq. (fO| shows that the quantum angle also controls the effect of population imbal- 
ance. The quantum formula, instead of including in the quasi-probabilities a factor: 



2^/NjV^ 

Na + Np 



(29) 



contains inside the integral an oscillating function cos [{Na — Np) A]. Figs. [5] and [5] illustrate the effects 
of population imbalance. 

To see what happens when the effect of A is cancelled, let us set A = in Eq. p^ . We write (l±cos A) 
in terms of sine and cosine of the half angle, expand these in exponentials, expand the binomials, and 
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integrate. The result is: 
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No oscillation then takes place; for instance, the case of = Np = 25 was already shown in Fig. [51 
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2.4 Pair-probability formulation 



When Na = Np — 1, the Hong-Ou-Mandel result is that the probabilities of having n particles in the 
detector 1 and and 2 — n in detector 2 is: 



(31) 



with n = 0, 1, 2. Now, a natural question is: can we consider that the distributions obtained above 
can be interpreted as those that one would be obtained by repeating the Hong-Ou-Mandel experiment a 
sufficient number of times, and accumulating the counts in each detector? 

First consider the case when Na = iV^ = N/2. We pick the first pair, and it produces either 2 
particles on the left or none. Then the second pair does the same. We continue until we have mi/2 pairs 
on the left and 1712/2 pairs on the right. We consider that, on the left we have mi/2 filled pair slots and 
1112/2 empty pair slots, which could have been selected in any order. We interchange the slots, while not 
counting interchanges of the empty slots among themselves and the filled among themselves. We have 
then (A^/2)!/[(mi/2)!(TO2/2)!] different ways for getting the (mi,m2) probability, which is: 



pipair) iN/2y. 

^ (mi/2)!(m2/2)! 



[(i + (-iro] 



(32) 



where the last factor ensures that there are an even number of particles on each side. For 50 particles 
the result is shown in Fig. [TUl 




Figure 10: The pair probability for iV„ = 25, Np = 25. 



We see that the distribution obtained in this way still has the expected odd-even behavior (actually 
almost by construction), but that it does not reproduce the characteristic 1 / ^mim2 shape for the enve- 
lope. Instead, it tends to concentrate the most likely results around mi ~ 7712 ~ N/2, which is natural: 
if we have independent scattering events of pairs into either channel, chosen randomly, one expects that 
the most likely values will be equal in both outputs. 
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By contrast, the curve of Fig.[5]has its maxima at mi — and nii = N, which is a completely different 
behavior, and indicates the effect of bosonic statistics (bunching into one channel). The distribution is 
no longer peaked at the center but spreads towards both sides. As Fig. [2] shows, this behavior can be 
explained if we add a new ingredient, a relative phase. But, since Fock sources do not have any initial 
phase, this can be understood as a result of a spontaneous choice of a relative phase by the two sources 
under the effect of quantum measurement }llj ; since this phase is completely unknown, an average over 
all possible values is taken. 

Now, if we have an excess of particles on one side, for example, = Njj , then we can assume each 
P particle is paired with an a particle and the M extras appear anywhere in the sequence of selections as 
singles. On the left side, we have Np pair slots with / of them filled and e of them empty. Of the single 
slots on the left, s are filled and o are open. Then we have 



Np = f + e 
Na- Nfi = s + o 
mi = 2/ + s 



(33) 



These can be solved to give: 



Nfi- - (mi - s) 



(34) 



Then among the pairs and singles on the left we can rearrange in / + e + s + o = Na total ways with 
rearrangements among the same kind of slots not counting to give a probability 
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..0 -Np- s)\s\ {Np - i(mi - s))\ 

The probability vanishes if mi — s is odd. The result for one extra particle is shown in Fig. 1111 The 
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Figure 11: The pair probability for iV„ = 26, = 25. 



results are then very different from those of Fig. [71 If there are two excess particles the steps on the side 
of the peak are completely smoothed out. The conclusion of that, in this case, the model of independent 
and repeated Hong-Ou-Mandel scatterings does not provide a good representation of the phenomenon at 
all. 
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3 Conclusion 



A single beam splitter cannot exhibit quantum non-local effects; violating local realism requires the 
combination of several such devices to form interferometers Nevertheless, here we have seen one 
beam splitter is sufficient to obtain interesting quantum effects, provided it receives Fock states at its two 
inputs; these effects are similar to the "population oscillations" predicted in more elaborate cases \TE[ I19| . 
The oscillations are related to the Hong-Ou-Mandel effect, but they cannot be understood as a simple 
juxtaposition of many separate two-photon-experiments. Actually, many-boson effects take place as a 
consequence of quantum statistics, which can be understood as a consequence of the tendency of two Fock 
states to acquire a relative phase under the effect of quantum measurement. Since this phase is completely 
unknown, the characteristic dependence shown in Fig. [5] results, onto which quantum oscillations are 
superimposed. Experimentally, the major difficulty for observing these effects is the production of Fock 
states with well-defined populations. Nevertheless, the experimental techniques that have been developed 
for Bose-Einstein condensates in ultra-cold gases seem well suited to planning experiments with input 
states that contain for instance a few tens of bosons. 

Note added in proofs: the authors have recently become aware of Ref. [20J, which gives a study of the 
localization of phase obtained by measurements of particles at the output of a single beam splitter. The 
theoretical treatment is similar to ours but this reference does not assume that all particles are detected; 
the oscillations we have discussed (generalization of the Hong-Ou-Mandel effect) do not appear. 
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